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We study supersymmetry breaking in the lattice N=l Wess-Zumino model by the world-line 
path integral algorithm. The ground state energy and supersymmetric Ward identities are exploited 
to support the expected symmetry breaking in finite volume. Non-Gaussian fluctuations of the 
topological charge are discussed and related to the infinite volume transition. 

PACS numbers: 11.30.Pb, ll.10.Ef, ll.10.Kk 



o 
o 

(N 



in 

(N 



(N 
O 

m 
o 
m 
O 

<D 
43 



X 



An important issue in the study of lattice supersym- 
metry is the analysis of the phase transition associated to 
spontaneous supersymmetry breaking [1]. The problem 
is particularly difficult in 1+1 dimensions where non per- 
turbative effects can be dominant [2]. In this Report, we 
address the determination of the phase diagram in the 
N=l Wess-Zumino (WZi) model. The transition must 
be analyzed by accurate sampling of the ground state 
and, to this purpose, there exist well established specific 
numerical techniques, e.g. Green Function Monte Carlo 
(GFMC) [3]. In a recent series of papers [4], the model 
WZi has been studied by GFMC with adaptive trial 
wave function optimization. The main source of system- 
atic error turns out to be the finite number K of quan- 
tum states that represent stochastically the ground state. 
The K — > oo extrapolation is delicate because the ana- 
lytical form of the asymptotic corrections is not known. 
For these reasons, it is worth exploring alternative meth- 
ods like the so-called world- line path integral ( WLPI) [5] , 
more common in the study of high energy lattice mod- 
els. In the framework of WLPI quantum expectation 
values are evaluated at finite temperature T = 1//3 and 
the ground state projection (/3 — > oo) is numerically esti- 
mated. In one spatial dimension, fermion configurations 
are sampled by local moves of the fermion lattice world 
lines. Previous WLPI studies of WZ X considered small 
lattices at moderately low temperatures [6] . A signal for 
the transition was given in terms of a small ground state 
energy compatible with zero within the achieved preci- 
sion. This interpretation is problematic since analytical 
arguments [7] as well as GFMC data [4] support SUSY 
breaking in finite volume for the considered model. A 
first aim of this Report is that of clarifying the results 
obtained with the WLPI method. We extend the analy- 
sis by including SUSY Ward identities and discussing the 
removal of systematic errors on larger lattices and lower 
temperatures. The topological charge is also measured 
and non-Gaussian fluctuations are found. Recently, the 
fluctuation spectra of global observables has been pro- 
posed as a monitor for criticality even relatively far away 
the transition [8]. Indeed, we shall motivate the con- 
jecture that the topological charge plays this role in the 
WZi model. 



The (1,1) SUSY algebra in two dimensions has two 
fcrmionic generators Qi,2 and three bosonic ones: the 
components of the two-momentum (P^P 1 ) and a cen- 
tral charge Z. The non trivial algebra is 



{Qa, Qb} = 2(P° 1 + PV + Za 3 ) ab , 



(1) 



where a 1 are Pauli matrices. The WZi model realizes the 
above algebra on a real chiral multiplet with a real scalar 
component tp and a Majorana fermion with components 
ipijj. The supercharges are 



Ql,2 = dx 



P ip 



1,2 



dip 
dx 



± v(<p) V 



(2) 



where p{x) is the momentum operator conjugate to 'fix) 
and V(ip) is an arbitrary function called prepotential in 
the following. The central charge is a topological quan- 
tum number [9] related to the existence of soliton sectors. 
On the lattice there are no continuous translation and we 
preserve a SUSY subalgebra [10]. We choose one of the 
supercharges, say Qi, build a discretized version Q and 
define the lattice Hamiltonian to be H = Q 2 . By Monte 
Carlo, we look for the eigenspace of Q 2 with lowest eigen- 
value. The simulation explores all the topological sectors: 
the Q-invariant states with Z = are expected to be lat- 
tice versions of fully supersymmetric ground states. In 
principle, there are also Q-invariant states with Z =/= 0, 
i.e. BPS states [9]. The transition between an empty 
Q 2 = eigenspace and a non trivial one is a weak (ac- 
tually halved) definition of SUSY breaking that we shall 
adopt in this Report [4]. 

The explicit lattice model is built by considering a spa- 
tial lattice with L sites and open boundary conditions. 
On each site we assign a real scalar field <p n with its 
conjugate momentum p n such that [p n ,(p m ] — —iS n ,m- 
The associated fermion is a Majorana fermion ip an with 

a = 1,2 and {lp a , n , 1pb,m} = ^a,b^n,m , V'a.n = V'a.n- Tne 

discretized fermionic charge and lattice (central) topo- 
logical charge are 



Q = zZ 



Pntp 



1,71 



Vn-l 



V(lfin) I ^2,n 



. (3) 
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n=l 



<Pn+l — <Pn-l 



V(<p n ) 



(4) 



For a smooth configuration the continuum limit of Z is a 
surface term. Following [6] we replace the two Majorana 
fermion operators with a single Dirac operator c satisfy- 
ing canonical anticommutation rules, i.e., {c n ,c m } = 0, 

{4i c m} = 0, {c n , C-ln} — Sn,m'- 



i>{l,2},n = 



tDl±i ( j ±iCn) 



The Hamiltonian is H = Q 2 = Hb + Hf with 



2 P " + 2 t 2 + V ^ n) 

n=l k 



L f 1 

#F = X! i -9 ( C n c «+1 + <i +1 C n ) + 
n=l L Z 



(5) 

(6) 
(7) 



c t c „-_ 



_ff conserves the total fermion number 

L 



N f = J2 Nn 



c n c m 



(8) 



n=l 



and can be examined separately in each sector with fixed 
Nf. The simplest observable that we consider is the lat- 
tice ground state energy Eo = minspec(<5 2 ) > 0. Ad- 
ditional quantities are related to global supersymmet- 
ric Ward identities as follows. If the vacuum \Q) obeys 
Q\Q) = 0, for any operator A we have (f2| {Q, A} \ = 0. 
In particular, taking A — J2n=i F(<Pn) ^2,n, we obtain 



w = (01 { n<p« 



n=l 



fn+1 — <Pn-l 



+ V{<p r , 



+ 



(9) 



+F'(^)(-l)«(4c„-l/2)}|0)=0. 



A set of independent Ward Identities is obtained by con- 
sidering F(ip) = ip k . We shall call Wk the corresponding 
quantities. In summary, Q\Q) = <-» Eo = and also 
<3|f2) = — > IFfc = 0. The world-line path integral com- 
putes the expectation value of an operator A at inverse 
temperature (3 as Tr (A e _/3H )/Tr e~ f3H . As usual the 
trace is computed by Trotter splitting of e~P H with time 
spacing e = 0/T 



-0H 



e — = (e- £UB e- £U ^e- £ ^) J +0{e z ) 
where H e J o) = ^„ even(odd) H F,n, 

H F , n = -\{cic n+1 + h.c.) + U-l) n V'( Vn ) ^N n 

Hn + l{-l) n+1 V'( Vn+1 )(N n+1 - 1 - 



(10) 



(11) 



with £„ = 1 on sites and L—l and 1/2 elsewhere. The 
terms Hp n with even n commute among themselves as 
well as the terms with odd n. Complete sets of inter- 
mediate states are inserted between the various factors 
in Eq. (10). The space-time lattice is shown in Fig. (1). 
The weight of a fermion configuration is computed by 
multiplying the weights associated to the shaded squared 
plaquettes. 



z 



z 



allowed move 



FIG. 1. Space time lattice. A world-line configuration with 
4 fcrrnions is shown. 



The finite temperature expectation values of the en- 
ergy, Z and the Ward identities can be readily measured. 
We choose to evaluate the kinetic bosonic energy by the 
Virial theorem. A full sweep is a heat bath on the un- 
shaded plaquettes with the elementary moves shown in 
Fig. (1) followed by a multi-hit Metropolis update on the 
field ip configuration. A detailed discussion of such algo- 
rithms can be found, e.g., in [5]. 

For the numerical analysis, we focus on the specific 
prepotential V(ip) = Xf 2 + c that despite its simplicity 
already shows all the interesting features of the symme- 
try breaking. At the classical level, SUSY is broken for 
c > 0. For negative c, the zeroes of V are associated to 
classical SUSY vacua and their fate after radiative cor- 
rections must be investigated. For this particular V, the 
analysis of the quantum model in the continuum can be 
pursued in full details and leads to the following con- 
clusions [7]. In finite volume, tunneling effects lift the 
ground state energy to a positive value for all (positive 
and negative) c. In infinite volume and sufficiently large 
c < 0, there can be inequivalent vacua with (ip) ^ 0. The 
fermion acquires a mass and the Goldstino required for 
spontaneous breaking is absent. The phase diagram of 
the lattice model is studied in the plane (A, c) of adimen- 
sional couplings. The continuum limit is obtained when 
A — > 0. The Rcnormalization Group trajectories with 
constant Physics are: 



c — ^-A log A = constant. 



(12) 



In the scaling region, SUSY is expected to be restored 
below a certain critical c. If the physical volume is not 
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large enough, the transition is smoothened and a resid- 
ual positive energy remains (even for quite negative c) 
as a consequence of the above mentioned SUSY break- 
ing in finite volume. The data analysis must be done 
after the removal of the systematic errors associated to 
e and (3. Expectation values are obtained only after the 
combined limit e — > followed by (3 — ► oo. This ex- 
trapolation makes the determination of the transition a 
rather expensive problem. In the first part of our anal- 
ysis, we consider the e and (3 dependence of the various 
observables. To this aim, we perform explicit simulations 
with L = 22 and the two values A = 0.35, 0.5 that are 
reasonably in the scaling region [4]. We take measures 
at T = 50, 100 with 2 < (3 < 8. The number of full 
updates is typically 10 6 . We apply an high number of 
Metropolis hits (~ 50) in order to match the accuracy of 
the fermionic heat bath sampling. We always look at the 
sector with Nf = L/2 where the ground state is expected 
to lie [4] . In Fig. (2) we show Eq at various e = (3 /T and 
(3; we also show the results of a quadratic extrapolation 
to e — > followed by an exponential (plus constant) one 
to estimate the (3 — > oo limit. 
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FIG. 2. Ground state energy per site Eo/L. The left graph 
shows the intermediate results at variable e and (3. The right 
graph shows the extrapolation to e — > followed by (3 — > oo. 




FIG. 3. Intermediate results for Z and Wo, 2, 4- Empty (full) 
symbols are obtained on the 22 x 50 (22 x 100) lattice. 

The systematic error due to e increases with (3 at fixed 
T simply because the corresponding e is bigger. The limit 
e — > is absolutely necessary and, in particular, restores 
positivity of Eq at large (3. As c is decreased, the ex- 
trapolated Eq reaches a positive plateaux. This residual 
value is slightly higher at A = 0.35, due to the larger 
finite volume effects. The beginning of the flat region is 
at c(0.5) ~ —0.4 and c(0.35) ~ —0.3, in agreement with 
Eq. (12). In Fig. (3) we show the intermediate data for 
Z and Wo, 2- It is interesting to note that Wo is quite 
independent on e at the explored values and is thus the 
most safe quantity to be analyzed. Fig. (4) collects the 
extrapolated data. 
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FIG. 4. Extrapolated results for Z, Wo,2,4 obtained from 
the data shown in Fig. (3). 

The central charge Z is definitely negative with a sharp 
variation around c ~ —0.3. If SUSY is to be recovered in 
the L — > 00 limit, we expect BPS states to be present to- 
gether with eventual Q = Z = states. The sign of Z is 
easily checked: in the WZi model we have Q\ = P° + Z 
and P° = Q\ + Q\ is a positive operator. Hence, states 
in the lowest Q\ eigenspace prefer Z < and states with 
Qi — have rigorously Z < 0. The Ward identities are 
all non zero although Wo and W2 approach smaller val- 
ues as c decreases. Again Wo flattens below c ~ —0.3. 
We conclude that SUSY is broken (at L = 22) for all c 
in agreement with the GFMC analysis of [4] and the an- 
alytical arguments [7]. The SUSY restoration observed 
in [6] must be explained in terms of the moderate sta- 
tistical errors. Indeed, as we have shown, the absolute 
values of E are quite small and arise from large can- 
cellations between contributions from Hb and Hp. The 
search for SUSY restoration requires a suitable finite size 
scaling analysis of data at rather larger L. The critical 
point is to be searched around c = — 0.4(— 0.3) for the 
considered A = 0.5(0.35). Indeed, the presented results 
suggest some kind of transition there. 
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In the final part of this Report we want to show how a 
consistent additional information is provided by the fluc- 
tuations of the central charge Z. In the continuum, Z is 
a global quantity depending on the boundary conditions. 
On the lattice, its geometrical meaning is not clear due 
to the quantum roughness of the configurations. An his- 
togram analysis of Z reveals that definite non-Gaussian 
fluctuations are present. The density of the normalized 
variable ( = (Z - (Z))/({Z 2 ) - (Z) 2 ) 1 / 2 turns out to be 
consistent with the generalized Gumbel one [11] 



P(C) = AAcxp[a (K + v - e^ +v )\ 



(13) 



where N, /U and v are functions of a fixed by normal- 
ization, zero mean and unit variance. Fig. (5) shows the 
best fit coefficient a(c) and two sample spectra. 



normal distribution. In the broken phase Z is expected to 
vanish, but an inspection of the MC histories shows that 
rare configurations appear with unphysical kinks respon- 
sible for the tail of P(C) at negative C < 0. Further anal- 
ysis is needed to determine the precise functional form 
of Z, particularly at large \(\, and to clarify its relation 
with the L = oo SUSY transition. However, Fig. (5, a) 
demonstrates clearly that the study of P(() is in prin- 
ciple rich of information. It is a pleasure to thank M. 
Campostrini and A. Feo for discussions about the WZi 
model and E. Alfinito for bringing the role of extremal 
statistics to our attention. Financial support from INFN, 
IS-RM42 is acknowledged. 
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FIG. 5. Analysis of Z (A = 0.35). The graph (a) shows 
the coefficient a (defined in the text) as a function of c. The 
graphs (b,c) show two sample fluctuation spectra together 
with the best fit with Eq. (13). 

In Fig. (5, a) a change of behaviour occurs around 
c* ~ —0.3. When c > c* , a(c) is small and weakly 
dependent on both c and (3. Below c* , a(c) increases 
rapidly as c is decreased and the slope is larger as (3 in- 
creases. The values of the fitted /i, v for large a are such 
that y/a[i — ► 1 and ^/av — > 0. For a — > oo, the lim- 
iting form of Eq. (13) converge to a normal Gaussian. 
Indeed the spectrum asymmetry decreases as c crosses 
c*. For instance, Fig. (5, c) is almost consistent with a 
simpler Gaussian fit. Similar features are discussed in [8] 
where the particular value a = n/2 is observed and re- 
lated to the joint effect of finite-size, strong correlation 
and self-similarity. In the WZi model, the quantum sym- 
metry breaking is not explained in terms of self-similar 
structures and a is just a measure of the deviation from 
Gaussian statistics. It is tempting to interpret Fig. (5) as 
another signal for a finite size critical c*(L) with SUSY 
restoration in infinite volume at c — c*(oo). At a quali- 
tative level, in the symmetric phase soliton states appear 
with a definite Z and the MC measures fluctuate with a 
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